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ABSTRACT 

We report on the Swift monitoring of the candidate supergiant fast X-ray transient (SFXT) 
IGR J 1 64 18 — 4532, for which both orbital and spin periods are known (~ 3.7 d and ~ 1250 s, 
respectively). Our observations, for a total of ~ 43 ks, span over three orbital periods and 
represent the most intense and complete sampling of the light curve of this source with a 
sensitive X-ray instrument. With this unique set of observations we can address the nature 
of this transient. By applying the clumpy wind model for blue supergiants to the observed 
X-ray light curve, and assuming a circular orbit, the X-ray emission from this source can be 
explained in terms of the accretion from a spherically symmetric clumpy wind, composed of 
clumps with different masses, ranging from ~ 5 x 10 16 g to 10 21 g. Our data suggest, based 
on the X-ray behaviour, that this is an intermediate SFXT. 
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1 INTRODUCTION 

The X-ray transient IGR J 1641 8 — 4532 was discovered by IN- 
TEGRAL during observations of the black hole X-ray transient 
4U 1630-47 on 2003 February 1-5 (Tomsick et al. 2004), with a 
20-40 keV flux of 3 x 10 -11 erg cm -2 s -1 . Sguera et al. (2006) re- 
ported on further INTEGRAL observations, taken on 2004 February 
26, during which fast (~ 1 hr) X-ray outbursts were observed that 
peaked at ~ 80 mCrab (20-30 keV). This behavior led Sguera et al. 
(2006) to propose that IGR J 1 64 18 — 4532 isacandidate supergiant 
fast X-ray transient (SFXT). SFXTs (e.g. Sguera et al. 2005) are 
high mass X-ray binaries (HMXBs) associated with OB supergiant 
stars via optical spectroscopy that show an X-ray dynamic range of 
3-5 orders of magnitude. Indeed, they are characterized by bright 
(peak luminosities of 10 36 -10 3 ' erg s -1 ), short (a few hours, as ob- 
served by INTEGRAL, Sguera et al. 2005; Negueruela et al. 2006) 
X-ray outbursts significantly shorter than those of typical Be/X- 
ray binaries, and a quiescent luminosity of ~ 10 32 erg s _1 (e.g. 
in’t Zand 2005; Bozzo et al. 2010). 

In their systematic analysis of the INTEGRAL observations 
from 2003 to 2009, Ducci et al. (2010) detected 23 outbursts, for 
an activity duty cycle of ~ 1 % (one of the highest among the 14 
SFXTs and SFXT candidates they examined) and fluxes ranging 


between 1.3 x 10 10 erg cm 2 s 1 and 4.8 x 10 10 erg cm 2 
s _1 (18-100 keV). 

XMM-Newton observations (Walter et al. 2006) showed a 
heavily absorbed (Vh = (l.OdtO.l) x 10 23 cm -2 ) pulsar (P sp m = 
1246±100 s) with an unabsorbed flux of 1.3 x 10 _1 ° erg cm -2 s _1 
(2-100 keV) in a fit performed in conjunction with the INTEGRAL 
data. 

Chaty et al. (2008) proposed, based on the XMM-Newton po- 
sition, 2MASS J 1 64 15078 — 4532253 as the best NIR counterpart 
candidate; they also fitted the spectral energy distribution (SED), 
obtaining a temperature of 32 800 K for the massive compan- 
ion, consistent with an OB spectral classification. SED fitting by 
Rahoui et al. (2008) yielded Av = 14.5 mag, and a 08.5 spectral 
type; depending on the luminosity class they derive a distance of 
4.9 (V), 8.3 (III), or 1 3 kpc (I). 

A periodicity of ~ 3.75 d was discovered by Corbet et al. 
(2006) in the Swift / Burst Alert Telescope (BAT, Barthelmy et al. 
2005) (3.753 ± 0.004 d, MJD 53360-53630) and RXTE1 ASM 
(3.7389 ± 0.0004 d, MJD 50091-53810) light curves, which was 
interpreted as the orbital period of a binary system. The hint of a to- 
tal eclipse at phase 0.5 observed in the BAT data is consistent with 
a supergiant companion. Levine et al. (2010) refined the orbital pe- 
riod value to 3.73886±0.00028 d, by using the RXTEI ASM data up 
to MJD 55224. The observed values of the orbital and pulse periods 
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Table 1 . Summary of the 5JW/I/XRT observations. 


ObsID 

Start time (UT) 

End time (UT) 

Exp. 

(s) 

r 

00031929001*' 

201 

1-02-18 

01:59:43 

201 

1-02-18 

03:43:26 

1951 

0.27 

00031929002 

201 

1-03-04 

00:03:19 

201 

1-03-04 

05:01:57 

2001 

0.01 

00043 105001 c 

201 

1-06-27 

05:14:54 

201 

1-06-27 

05:23:56 

537 

0.80 

00031929003 

201 

1-07-13 

09:46:55 

201 

1-07-13 

14:43:57 

4756 

0.16 

00031929004 

201 

1-07-14 

12:39:39 

201 

1-07-14 

19:31:56 

4832 

0.47 

00031929005 d 

201 

1-07-15 

06:42:44 

201 

1-07-15 

23:00:58 

273 

0.72 

00031929005 

201 

1-07-15 

06:42:50 

201 

1-07-15 

23:01:57 

1078 

0.72 

00031929006 d 

201 

1-07-16 

10:15:02 

201 

1-07-16 

23:05:36 

458 

0.01 

00031929006 

201 

1-07-16 

10:15:05 

201 

1-07-16 

23:06:58 

379 

0.01 

00032037001 

201 

1-07-18 

21:28:28 

201 

1-07-19 

00:51:31 

1947 

0.62 

00031929007 

201 

1-07-19 

01:52:41 

201 

1-07-19 

23:19:57 

2673 

0.77 

00032037002 

201 

1-07-20 

05:18:42 

201 

1-07-20 

15:22:58 

3483 

0.01 

00032037003 

201 

1-07-21 

10:06:18 

201 

1-07-21 

15:27:54 

3761 

0.31 

00032037004 

201 

1-07-22 

15:19:43 

201 

1-07-22 

18:44:57 

2812 

0.62 

00032037005 

201 

1-07-23 

11:48:19 

201 

1-07-23 

15:23:56 

4127 

0.85 

00032037006 

201 

1-07-29 

12:13:30 

201 

1-07-29 

23:34:57 

4650 

0.50 

00032037007 

201 

1-07-30 

01:01:30 

201 

1-07-30 

09:27:56 

4287 

0.63 


a Mean phase referred to T epoc h = MJD 53560.2 (Levine et al. 2010). 
b Corresponding to MJD 55610.08. 
c Serendipitous observation. 
d WT mode data. 

indeed place IGR J 1 64 18 — 4532 in a region of the Corbet diagram 
(Corbet 1986) where wind-accreting HMXBs generally lie. 

In this paper we analyze all the Swift data collected on 
IGR J 1641 8— 4532. The first set of observations comprises two 2ks 
target of opportunity (ToO) observations performed on 201 1 Febru- 
ary 18 and March 4, following a MAXI alert of a possible detec- 
tion of IGR .116418— 4532 on 2011 February 17 (preliminary re- 
sults were described in Romano et al. 201 lb). The second set is 
drawn from a monitoring campaign in 2011 July that spans over 
three orbital periods. This Swift monitoring represents the most in- 
tense and complete sampling along the orbital period of the light 
curve of this source with a sensitive X-ray instrument. 


2 OBSERVATIONS AND DATA REDUCTION 

Table 1 reports the log of the Swift/X - ray Telescope (XRT, 
Burrows et al. 2005) observations used in this paper. The Swift ob- 
servations of IGR J 1641 8 — 4532 during 2011 February and March 
(00031929001 and 00031929002) were obtained (Romano et al. 
2011b) as two ToO observations following a MAXI alert 
(5609524433, later retracted, as due to GX 340+0) of a possible 
detection ofIGRJ16418— 4532 on 2011 February 17. The data col- 
lected in 201 1 July were obtained as a ToO monitoring program of 
9 scheduled daily observations, each 5ks long, starting on 2011 
July 13. This strategy was devised in order to cover the light curve 
in different orbital phases while maintaining the observing time per 
day reasonably short, thus not hindering gamma-ray burst (GRB) 
observations. However, most observations were cut short by the 
occurrence of GRBs and the discovery of a new soft gamma-ray 
repeater (starting from obsID 00031929005, see Table 1). Obser- 
vation 00032037001 was performed to determine the best offset 
from the intended target in order to alleviate the problem of single- 
reflection rings from the nearby source GX 340+0 (these rings 


marginally affect the source while in the faint states, so that the 
analysis of the earlier segments was still feasible, with an appro- 
priate choice of extraction regions) A new pointing (hence a new 
obsID) was adopted afterwards. In our work we also include ob- 
servation 00043105001 that serendipitously covers the region of 
IGR J 1641 8 — 4532. The 2011 July campaign lasted 18 days di- 
vided in 12 observations for a total on-source exposure of ~ 39 
ks. 

For our binary orbital period search and analysis we retrieved 
the BAT ‘orbit-by-orbit’ light curves (averaged over Swift ’ s orbital 
period of ~ 90 m, 15-50 keV) covering the data range from 2005 
February 12 to 2011 July 12 (MJD range 53413-55754) from the 
BAT Transient Monitor (Krimmetal. 2006, 2008) page 1 . These 
data were further screened to exclude bad quality points (quality 
flag 1 and 2) and referred to the solar system barycentre (SSB) by 
using the task EARTH2SUN. 

The XRT data were processed with standard procedures 
(XRTPIPELINE vO.12.6), filtering and screening criteria by using 
FTOOLS. Both WT and PC events were considered. The selec- 
tion of event grades was 0-2 and 0-12, for WT and PC data, re- 
spectively (Burrows et al. 2005). Source events were accumulated 
within a annular/circular region (depending on whether pile-up cor- 
rection was required or not, respectively) with an outer radius of 20 
pixels (1 pixel ~ 2.36”); and background events were accumulated 
from source- free regions far from the contamination of the single- 
reflection rings from GX 340.0+0. For our timing analysis we con- 
verted the event arrival times to the Solar system barycentre with 
the task BARYCORR. Light curves were created for several values 
of signal-to-noise ratio (SNR) and number of counts per bin; all 
were corrected for PSF losses, vignetting and background. For our 
spectral analysis, we extracted events in the same regions as those 
adopted for the light curve creation; ancillary response files were 
generated with XRTMKARF, to account for different extraction re- 
gions, vignetting, and PSF corrections. We used the latest spectral 
redistribution matrices in CALDB (201 10705). 

The UV/Optical Telescope (UVOT, Roming et al. 2005) ob- 
served IGR J 1641 8 — 4532 simultaneously with the XRT with the 
‘Filter of the Day’ (FoD), i.e. the filter chosen for all observations 
to be carried out during a specific day in order to minimize the fil- 
ter wheel usage. The exception is observation 00032037007, which 
was requested in the u filter, as opposed to FoD, to attempt to ob- 
tain a detection in that band. The data analysis was performed using 
the uvotimsum and UVOTSOURCE tasks included in the FTOOLS. 
The latter task calculates the magnitude through aperture photom- 
etry within a circular region and applies specific corrections due to 
the detector characteristics. 

All quoted uncertainties are given at 90% confidence level 
for one interesting parameter unless otherwise stated. The spectral 
indices are parameterized as F u oc u~ a , where F v (erg cm” 2 s” 1 
Hz” 1 ) is the flux density as a function of frequency v\ we adopt 
T = a + 1 as the photon index, N(E) oc E~ r (ph cm” 2 s” 1 
keV” 1 ). 

3 RESULTS 

3.1 Orbital period from BAT data 

For the orbital modulation, we considered the BAT data in the 
time range MJD 53414-55754. With a standard folding analy- 

1 http://swift.gsfc.nasa.gov/docs/swift/results/transients/ 
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Figure 1. BAT light curve folded at P or b = 3.73886 d and T ep och = 
MJD 53560.20000. 

sis of the Solar system barycentred light curves, we measured 
an orbital period of 3.740 ± 0.002 days (1 -a error), which is 
fully consistent with the value obtained by Levine et al. (2010), 
Porb = 3.73886 ± 0.00028 d. Therefore, hereafter we shall adopt 
the more precise measure by Levine et al. (2010), and we shall refer 
to Tepoch = MJD 53560.20000 (Corbet et al. 2006). Fig. 1 shows 
the BAT light curve folded over this period. 

The BAT data show an eclipse centred at ij> ~ 0.55. The 
depth is consistent with a total eclipse, and the eclipse FWHM 
duration is 0.17 ± 0.05 of the orbital period. The folded light 
curve is quite flat; the peak-to-trough pulsed fraction, PF p t = 
(T m ax T n i n ) / ( /''ncLx ± F m in), where Tmax and Tmin are the 

observed background-subtracted count rates at the peak and at the 
minimum, is 78 ± 12 % (1-cr error). 

3.2 X-ray position 

We used 20.4ks of PC mode data and simultaneous SwiftlUVCtT 
images to obtain an astrometrically-corrected position (see 
Evans et al. 2009; Goad et al. 2007) of: RA(J2000) = 250.46104, 
Dec(J2000) = —45.54092, which is equivalent to: RA(J2000) = 
16 h 41 m 50165, Dec(J2000) = -45° 32' 27'.' 3, with an uncer- 
tainty of l'/9 (90% confidence level) 2 . This position is 2'/5 from 
the source 2MASS J16415078-4532253 and 4"3 from the XMM- 
Newton position (Walter et al. 2006), and it allows us to confirm the 
2MASS source as the optical counterpart of IGR .1 16418—4532. 

3.3 UVOT optical counterpart 

At the XRT position of IGR J 1641 8 — 4532, no detection was 
achieved in any UV band down to a limit of wl > 21.78, 
m2 > 21.66, w2 > 21.89 mag, respectively. A marginal de- 
tection (2.9cr) was obtained in the u band, at 22.28 ± 0.40 mag. 
These magnitudes are on the UVOT photometric system described 
in Poole et al. (2008), and are not corrected for Galactic or intrinsic 
extinction. At this position, the column due to Galactic absorption 
is 1.59 x 10 22 cm -2 (Kalberla et al. 2005). By adopting the con- 
version into Ay (Predehl & Schmitt 1995), Nh = 1.79 x 10 21 Ay 

2 This new UVOT-enhanced XRT position takes advantage of the new cal- 
ibration of the XRT to UVOT detector co-ordinates of 20 1 1 August (Evans 
2011); see, also: http://www.swift.ac.uk/reprocessed.php . 



Time (s) 

Figure 2. Swifl/XKT 0.3-10 keV light curve of IGR J16418— 4532 dur- 
ing observation 00031929001, background-subtracted, corrected for pile- 
up, PSF losses, and vignetting (SNR = 5). Time is referred to the start of 
the observation. 

cm -2 , this corresponds to (Cardelli et al. 1989) A u = 14.2 mag in 
the UVOT u band. 


3.4 Spin period search 

We looked for evidence of the 1246 ± 100 s spin periodicity re- 
ported by Walter et al. (2006) in the Solar system barycentred XRT 
light curves. We only considered observations where the source 
was significantly detected (thus excluding obsIDs 00032037001 
and 00032037004). 

We searched each single observation longer than one spin pe- 
riod and the total combined dataset within ±300 s (corresponding 
to a ±3cr error reported by Walter et al. 2006) by applying both a 
Fourier transform and a Z 2 analysis (a method that does not re- 
quire the binning of the data; Buccheri et al. 1983). No statistically 
significant signal was detected either in the individual observations 
or in the combined dataset. With respect to the latter, the upper 
limit on the pulsed fraction, (defined as the semi-amplitude of the 
sinusoidal modulation divided by the mean count rate), computed 
according to Vaughan et al. (1994), is 42 % at the 99 % confidence 
level and for periods between 946 and 1546 s. 

This upper limit is only marginally compatible with the (64 ± 
10) % value measured by Walter et al. (2006). We note however 
that they were able to detect pulsations only during a fraction of 
their XMM-Newton data and in particular while the source was not 
flaring, and that our data were collected in a variety of states of 
brightness and across almost all orbital phases. So our upper limit 
and their measure may not be straightforward to compare. 

3.5 Light curves 

Fig. 2 shows the 0.3-10keV light curve of IGR .116418—4532 dur- 
ing the 2001 February 18 observation, binned in order to achieve 
a SNR of 5. The source reached ~ 5 counts s^ 1 during the ini- 
tial flare and then decayed through a series of smaller flares. The 
lowest point is at 0.16 counts s^ 1 . When binned at a SNR of 3 the 
dynamical range, during this observation, is of order 50. 

Fig. 3 shows the 0.3-10 keV light curve of IGR .116418—4532 
of the whole 2011 July campaign after bin-by-bin background- 
subtraction and corrections for pile-up, PSF losses, and vignetting, 
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Figure 3. Swift/XKT 0.3-10 keV light curve of IGR J 1 64 18 — 4532 during the 201 1 July monitoring program, background-subtracted and corrected for pile-up, 
PSF losses, and vignetting (SNR= 3). Downward-pointing arrows are 3 a upper limits. Different colours (see the electronic edition) mark different observations 
(see Table 1). The top axis reports the phase (with P Q rb — 3.73886 d and T epoc h = MJD 53560.20000). 
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Figure 4. Swift/XKT 0.3-10 keV flux light curve of IGR J1641 8-4532, folded at P orb = 3.73886 d and T epoch = MJD 53560.20000. The data were 
collected in 201 1 February and June (dark grey filled circles), and 2011 July (black filled circles). Downward-pointing arrows are 3cr upper limits. The open 
blue circle is a detection obtained by combining three observations (00032037001, 00032037004, and snapshots 1-3 of 00032037007). Also included are the 
2004 August XMM data (light grey curve). 


binned at SNR = 3. The light curve starts at phase 0.16 (assum- 
ing a period of P orb = 3.73886 d and an initial epoch T ep och = 
MJD 53560.20000) and monitors the source state through over 
three full periods. Superimposed on the long-term orbital modu- 
lation, which follows that seen in the BAT data, flaring is observed 
on short time scales, as shown in Fig. 2. This behaviour has been 
observed in most SFXTs (Sidoli et al. 2008; Romano et al. 2009b, 
2011a, 2010). 

The lowest point (as observed by XRT) was collected dur- 
ing the eclipse, by combining three observations (00032037001, 


00032037004, and snapshots 3 1-3 of 00032037007, see open blue 
point in Fig. 4) at 2.5 x 10~ 3 counts s -1 . This corresponds to an un- 
absorbed 0.3-10 keV flux of ~ 6.0 x 10 -13 erg cm -2 s _1 , and to a 
luminosity of 1.2xl0 34 erg s^ 1 (assuming the optical counterpart 
distance of 13kpc). The lowest point during the campaign outside 
the eclipse was recorded on MJD 55763.43 at 0.022 counts s _1 . 


3 Part of an observation with a continuous viewing of the target, with typi- 
cal timescales of 10 3 s. 
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Figure 5. Swift/XKT hardness ratio (red open circles, 30 counts per bin) 
as a function of count rate (top) and phase (bottom). The black points are 
weighted means. 

corresponding to an unabsorbed 0.3-10 keV flux of ~ 4.7 x 10” 12 
erg cm” 2 s” 1 , and to a luminosity of 9.5 xlO 34 erg s” 1 . As the 
peak count rate is reached on MJD 55610.08 at > 8 counts s” 1 
the observed dynamical range of this source is at least 370, (1400 
considering the points within the eclipse). 

Fig. 4 shows the 0.3-10 keV light curve, folded using P or b = 
3.73886 d and T epoc h = MJD 53560.20000. We adopted count rate 
to flux conversion factors derived from spectral fits of each obser- 
vation (Section 3.6) and subsequent interpolation across observa- 
tions where fitting was not feasible. We also include the ~ 10 ks 
2004 August 19 XMM-Newtonl EPIC-pn observation (Walter et al. 
2006), binned at 256 s, and converted into flux by adopting a 
single conversion factor obtained from fitting the mean spectrum 
with an absorbed power-law model (F = 1.40 ± 0.05, Nn = 
(10.5 ± 0.5) x 10 22 cm -2 , 1 -a errors; + =0.99/169 dof). The 
eclipse lasts A (f> ~ 0.2 and is centred at <f> ~ 0.55, as observed in 
the BAT data (Fig. 1). 

Fig. 5 shows the hardness ratio (HR defined as H/S, where H 
is the count rate in the 4-10 keV band and S is the count rate in the 
0.3 — 4 keV band), as a function of count rate (top), and orbital phase 
(bottom), where the light curve was binned to have 30 counts per 
bin (red open circles). A fit to a constant value yields a hardness 
ratio value of 1.01 ± 0.03 and xt = 0.85 for 130 degrees of free- 
dom (dof), so that the HR is shown to be roughly constant over two 
orders of magnitude in count rate. Furthermore, the HR were also 


Table 2. XRT spectroscopy of IGR J16418— 4532. 


ObsID 

A h 

r 

Flux“ 

L b 

X 2 /d°f 


( 10 22 cm” 2 ) 




C-stat (dof) c 

00031929001 

ft 77+2.32 
u. 1 f _ 1<7g 

1 qq+ 0-61 
A ' d **-0.54 


43 

0.532/18 

00031929002 

4 17+1-70 

1 Q4 + 0,51 
1 ,U 4 t -0.45 

11.58+; 39 

23 

0.845/17 

00043105001 

k 28”*" 3 ' 47 

O.ZS _2 54 

1 qc;+0-99 
-L.oa_o 87 

3-85^; 19 

8 

70.05(76) 

00031929003 

5.20+J ; 99 

i qc+0.27 
l.OO—o 24 

q Q7+0.65 
— 0.58 

17 

0.935/58 

00031929004 

5-16^2 81 

0 09” 1 ” 0 ' 87 

u.uy _ 0 75 

1 qk+ 0-21 

l.UD—o 17 

2 

169.46(141) 

00031929007 

3 93 + 1 - 24 
0.97 

1 04” 1 ” 0 ' 40 

7 62”*” 0 ' 72 
‘ - OZ _0.69 

15 

0.647/28 

00032037002 

4 59 "*" 1,23 

1 05+ 0,38 
l.UO-o 34 

ft oq+0.56 
-0.54 

13 

0.831/31 

00032037003 

5.65+; 3 ° 

0 62”*” 0 ' 62 
U-O^_o.55 

9 ‘aI +0-32 
z ' ol -0.31 

5 

1.114/16 

00032037005 

5-i4±i;S 

1 01 +0 - 45 
1,u± -0.40 

o oq+0.40 
9.8o_o 37 

8 

0.838/27 


a Unabsorbed 2-10 keV fluxes (10” n erg cm -2 s” 1 ). 
b 2-10 keV luminosities in units of 10 35 erg s — F at 13 kpc. 
c Cash statistics (C-stat). 


weighted across several phase intervals (black points); we cannot 
find evidence of variations of the HR as a function of phase. 

3.6 Spectral analysis 

We first extracted PC spectra for each segment in which a detection 
was obtained and a minimum of ~ 100 source counts were avail- 
able. We did not consider the WT data for our spectroscopic study 
as they were statistically poor and were heavily contaminated by 
the background from the single-reflection rings from GX 340.0+0. 
The data were rebinned with a minimum of 20 counts per energy 
bin to allow \ 2 fitting, with the exception of low-count-statistics 
spectra, in which Cash (Cash 1979) statistics and spectra binned 
to 1 count per bin were used, instead. The spectra were fit in the 
0.3-10 keV energy range with a single absorbed power law model. 
More complex models were not required by the data. The results 
are reported in Table 2. We note that the spectral parameters do not 
show significant variability within the large uncertainties, despite 
the large observed variations in flux throughout the campaign. 

Observation 00031929001 (shown as the dark grey filled cir- 
cles at phase 0.27 in Fig. 4, and detailed in Fig. 2), in particular, 
deserves special attention. Although it was performed as a follow- 
up of a false trigger (later attributed to the nearby GX 340+0), it 
nonetheless found the source in a very active state, the brightest of 
the whole campaign, with its average 2-10 keV unabsorbed flux of 
~ 2x 10” 10 erg cm” 2 s” 1 , corresponding to ~ 4x 10 36 ergs” 1 at 
13 kpc). The peak flux was ~ 2x 10” 9 erg cm” 2 s” 1 (~4x 10 37 
ergs” 1 ). An absorbed power-law model yielded a high absorbing 
column (Vh = (6.8l 2 'g) x 10 22 cm” 2 (in excess of the Galac- 
tic one, 1.59 x 10 22 cm” 2 ) and a photon index F = 1.4+® (see 
Table 2). 


4 DISCUSSION 

IGR J 1641 8 — 4532 is a candidate SFXT, based on the hard X-ray 
behaviour (Sguera et al. 2006), for which we know both the or- 
bital (P or bit ~ 3.7 d, Corbet et al. 2006) and the spin (~ 1250 s, 
Walter et al. 2006) periods. These values place IGR J 16418 — 4532 
in the region of the Corbet diagram (Corbet 1986) of the wind- 
accreting supergiant high-mass X-ray binaries. In this Paper we 
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have taken advantage of the short orbital period which makes 
IGR J 16418 — 4532 a very good source to monitor with the XRT, 
to study the X-ray properties of this source as a function of the 
orbital phase, in particular, to address the issue of its SFXT nature. 

The Swj/l/BAT and XRT folded light curves of 
IGR J 1641 8 — 4532 (Fig. 1, and Fig. 4) show the presence of 
a dip that is consistent with zero intensity. When interpreted as 
an eclipse, it can be used to infer the nature of the primary of 
the binary system. The BAT lightcurve of IGR J16418 — 4532 
(excluding the eclipse) does not clearly show a flux modulation 
with the orbital phase, suggesting that the eccentricity is not 
large. This is in agreement with the observation that in many 
X-ray binary systems with nearly circular orbits and short orbital 
periods (^ 4 days), the eccentric orbit produced by the explosion 
of the supernova is circularized on timescale of 10 6 yr (consistent 
with the evolution timescale of an O star), due to viscous tidal 
interactions (see e.g. Press et al. 1975). Thus, for simplicity, we 
assume a circular orbit (P or b = 3.73886 d). Let us adopt the 
typical values of stellar masses and radii (Martins et al. 2005) 
for a primary of spectral type 08.5 and luminosity classes I, III, 
or V, and for the neutron star a mass of ATns = 1.4 Mq, and a 
radius of Rns = 10 km. By combining Kepler’s third law with the 
relationship given by Rappaport & Joss (1983) for the radius of the 
primary R, 

(R/a) 2 = cos 2 i + sin 2 i sin 2 0 e , (1) 

where a is the semi-major axis of the orbit, i is the inclination, 
and 9 e is the eclipse semi-angle, we find that the observed Acp = 
0.2 is inconsistent with a luminosity class III or V primary. We 
therefore conclude that indeed the primary is an 08.51 supergiant 
star (with mass, radius, luminosity, effective temperature given 
by Martins et al. (2005), Mob = 31.5Mq, Rob = 21.4Rq, 
log Lob /Lq = 5.65, T e g = 32274 K). Therefore, the distance 
to IGR .116418-4532 is - 13kpc (Rahoui et al. 2008). 

During the eclipse of IGR J 1641 8— 4532, we detect a non- zero 
flux with the XRT (Lx ~ 3 x 10 34 ergs -1 , see Fig. 4). Detec- 
tions of HMXBs during the eclipse are quite common (see e.g. 
van der Meer et al. 2005; Clark et al. 1988; Sako et al. 1999) and 
can be ascribed to X-ray photons scattered by the stellar wind of 
the OB star and discrete fluorescent and recombination lines from 
highly ionized gas around the compact object. The decrease in X- 
ray flux, when the compact object is approaching the eclipse, is 
related to the increase of the column density along the line of sight, 
through the stellar wind produced by the OB star. From obser- 
vations of other wind-fed X-ray pulsars, such as 4U 1700—377, 
Vela X-l, OAO 1657—415 (see e.g. van der Meer et al. 2005; 
Sako et al. 1999; Audley et al. 2006) the absorption column density 
during the eclipse ranges from « 10 23 cm -2 up to ~ 10 24 cm -2 . 

We can compare the X-ray luminosity and absorption column 
density of IGR J 1 6418 — 4532 during the eclipse with those cal- 
culated by van der Meer et al. (2005) for 4U 1700—377, a wind- 
fed HMXBs with properties similar to IGR J 1 64 18 — 4532 (orbital 
period P or b = 3.4 d, companion star of spectral type 06.5Iaf). 
van der Meer et al. (2005) observed 4U 1700—377 with XMM- 
Newton during an eclipse; they found an X-ray luminosity of 
ss 7 x 10 34 ergs _1 (in agreement with the X-ray luminosity of 
IGR J 1641 8 — 4532 during the eclipse) and a column density of 
« 10 24 cm -2 . 

The XRT light curve (Figs. 2 and 3) shows a large dynamic 
range, up to 1-2 orders of magnitude within a single snapshot and 
about 3 orders of magnitude overall (even when the low-intensity 
state is excluded). While this dynamic range is lower than that ob- 



Figure 6. Observed cumulative luminosity distribution (solid black line) 
and relative uncertainties (dashed black lines), compared with a possible 
solution (red dot-dashed red line; see text). 


Table 3. Wind parameters values for IGR J16418— 4532. 


Parameter 

min. value 

max. value 

-A^tot 

4 x 10 -7 M 0 yr“ 1 

8 X 1O“ 7 M 0 yr- 1 

/ 

0.85 

0.95 

Voo 

800kms _1 

1300 km s" 1 

P 

0.8 

1.2 

C 

0.8 

1.6 

7 

-1.5 

-0.5 

M a 

5 x 10 16 g 

5 x 10 19 g 

M b 

5 x 10 19 g 

10 21 g 


served in the classical (Chaty et al. 2010) SFXTs, which reach 4-5 
orders of magnitude in dynamic range (e.g. IGR J17544— 2619), it 
is indeed typical of intermediate SFXTs (e.g., IGR J18483— 031 1). 
The 2011 February flare (Fig. 2 and dark grey filled circles at 
phase 0.27 in Fig. 4), in particular, reached ~ 2 x 10 -9 erg 
cm -2 s _1 (2-10 keV, unabsorbed), corresponding to a luminos- 
ity of ~ 4 x 10 3 ' ergs -1 , which is typical of SFXT outbursts as 
we have observed them with Swift/X RT (e.g., Romano et al. 2008, 
2009a; Sidoli et al. 2009b, a). 

On the other hand, while the typical SFXT light curve is 
characterized by a faint flaring activity with luminosity of 10 33 — 
10 34 erg s -1 (Sidoli et al. 2008; Romano et al. 2009b, 2011a), spo- 
radically interrupted by flares with duration of a few hours and 
peak luminosity of 10 36 — 10 37 ergs _1 , IGR J 1 64 18 — 4532 ex- 
hibits a flaring activity ranging from ~ 10 35 ergs _1 up to ~ 
4x 10 3 ' erg s _1 (also see Walter et al. 2006 and the light grey curve 
in Fig. 4). 

To reproduce the observed X-ray variability, we applied the 
clumpy wind model of Ducci et al. (2009), which we briefly sum- 
marize here. This model assumes that the wind of OB supergiants is 
inhomogeneous, composed of dense clumps surrounded by a hotter 
and homogeneous wind. Clumps follow a power-law mass distribu- 
tion 

= <2) 

where k = (1 — f) — M 1-1 ”) is the normaliza- 

tion constant, M c \ is the mass of a clump, and M a and Mb de- 
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fine the extrema of the mass range. Assuming a spherical geometry 
for the clumps, the initial clump dimension distribution is given by 
Nm c i oc clumps s -1 , where R c l is the radius of the clump. 
The clump mass loss rate M c i is related to the total mass loss rate 
Mtot by means of f — M c \/Mtot, the fraction of wind mass in the 
form of clumps. Clumps and the homogeneous wind follow the 13- 
velocity law v(r ) = Voo (l — 0.9983 ^ (Castor et al. 1975), 
where [3 is a constant in the range ~ 0.5 — 1.5, Rob is the radius 
of the supergiant and Voo is the terminal wind velocity. 

IGR J 1 64 18 — 4532 was intensely monitored by Swift/XRT for 
two orbital cycles (MJD 55755-55766). Each snapshot has a du- 
ration of ~ Iks, for a total of up to ~ 5ks per day, which is 
the typical timescale for the flare durations observed in SFXTs 
(Sguera et al. 2006). For this reason, we cannot establish the num- 
ber of observed flares and their duration, therefore we cannot use 
this information to compare the observed lightcurve with that cal- 
culated with the clumpy wind model of Ducci et al. (2009), as has 
been done in previous work (Ducci et al. 2009; Romano et al. 2010; 
Ducci et al. 2010). However, it is possible to compare the observed 
cumulative luminosity distribution out of eclipse with that calcu- 
lated with the clumpy wind model. For the former, the Swift/XRT 
count rate, measured in the 0.3-10 keV, was converted to the 0.1- 
100 keV luminosity using the average spectral parameters obtained 
by Walter et al. (2006) and Ducci et al. (2010), and assuming a dis- 
tance of 13kpc. For the latter, we adopted the physical parame- 
ters for the 08.51 primary and the neutron star secondary reported 
above, and the force multiplier parameters (see Ducci et al. 2009) 
obtained by Shimada et al. (1994) for a 08.51 star ( k = 0.375, 
a = 0.522, 5 = 0.099). 

Fig. 6 shows the Swift /X RT cumulative luminosity distribu- 
tion (black solid line) with its 90% uncertainties (black dashed 
lines). Fig. 6 also shows the cumulative luminosity distribution 
(red dot-dashed line) calculated with the clumpy wind model by 
adopting the following parameter values: mass loss rate Mtot = 
6 x 10 — ' Mq yr _1 , terminal velocity Voo = 1000kms _1 ,/3 = 1, 
fraction of mass lost in clumps / = 0.9, mass distribution power- 
law index £ si* 1.1, power-law index of the initial clump dimension 
distribution 7 = —1, minimum clump mass M a = 5 x 10 18 g and 
maximum clump mass Mb = 10 20 g. Further acceptable solutions 
can be found assuming the wind parameters in the ranges reported 
in Table 3. 

The wind parameters obtained for IGR .116418— 4532 are 
roughly in agreement with those of HMXBs previously studied (see 
Ducci et al. 2009; Romano et al. 2010; Ducci et al. 2010), with the 
exception of / and Voo ■ The obtained low value of the terminal ve- 
locity can be explained as follows: when the wind is highly ionized 
by the X-ray photons emitted by the neutron star, the high ioniza- 
tion modifies the dynamics of the line-driven stellar wind of OB 
supergiants (especially in close binary systems), leading to a re- 
duction of the wind velocity in the direction of the neutron star (see 
e.g. Stevens & Kallman 1990 and references therein). The mass 
loss rate is lower than the typical mass loss rate of 08.51 stars, 
which is of the order of « 2 x 10 -6 M 0 yr -1 (Vink et al. 2000), 
but in agreement with the hypothesis that the mass loss rates de- 
rived from homogeneous-wind model measurements with the Ha 
method are overestimated by a factor of 2-10 if the wind is clumpy 
(see e.g. Lepine & Moffat 2008). However, we cannot exclude the 
possibility that the X-ray variability of IGR J 164 18— 4532 is not to- 
tally due to the accretion of an inhomogeneous wind. In fact, other 
mechanisms could be at work, like e.g. transient accretion discs 
(Taametal. 1988) and intermittent accretion flow onto a neutron 
star (Lamb et al. 1977; Ducci et al. 2010). 


Our results, based on our intense monitoring of 
IGR J 1 64 18 — 4532, show that its X-ray properties (dynami- 
cal range, flaring activity and peak luminosity) are consistent with 
an intermediate SFXT nature. 
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